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ABSTRACT 

We identify near-infrared Kg band counterparts to i7ersc/ie/- ATLAS sub-mm sources, using a pre- 
liminary object catalogue from the VISTA VIKING survey. The sub-mm sources are selected from 
the H-ATLAS Phase 1 catalogue of the GAMA 9h field, which includes all objects detected at 250, 
350 or 500 /im with the SPIRE instrument. We apply and discuss a likelihood ratio (LR) method for 
VIKING candidates within a search radius of 10^^ of the 22,000 SPIRE sources with a 5cr detection 
at 250 iim.. We estimate the fraction of SPIRE sources with a counterpart above the magnitude limit 
of the VIKING survey to be Qq ^ 0.73. We find that 11,294 (51%) of the SPIRE sources have a best 
VIKING counterpart with a reliability R > 0.8, and the false identification rate of these is estimated 
to be 4.2%. We expect to miss ~ 5% of true VIKING counterparts. There is evidence from Z — J 
and J — Kg colours that the reliable counterparts to SPIRE galaxies are marginally redder than the 
field population. We obtain photometric redshifts for ~ 68% of all (non-stellar) VIKING candidates 
with a median redshift of z = 0.405. We have spectroscopic redshifts for 3147 (~ 28%) of the reliable 
counterparts from existing redshift surveys. Comparing to the results of the optical identifications 
supplied with the Phase I catalogue, we find that the use of medium-deep near-infrared data improves 
the identification rate of reliable counterparts from 36% to 51%. 
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1 INTRODUCTION 

The extra-galactic universe has been well observed in the op- 
tical and, to an extent, the near-infrared wavelength ranges 
for several decades; in contrast, only relatively recently 
have we started to carry out unbiased surveys in the sub- 
mm wavelengjths (ISmail. Ivison Blainlll997l : iMortier et alJ 
I2OO5I : iDevlin et alJ l2Q09l ) and we still have a limited un- 
derstanding of the sources responsible for the bulk of sub- 
mm emissions. Finding optical/near-infrared counterparts 
enables us to complement our knowledge with optical/near- 
infrared colours, photometric redshifts, and potentially fol- 
lowup with multi-object spectroscopy. 

Using near-infrared instead of optical wavelengths to iden- 
tify sub-mm sources allows us to probe galaxies out to 
higher redshifts {z > 0.5) where rest-frame visible bands 
are shifted to the observed near- infrared. This is especially 
important for the dusty galaxies expected to be detected 
by ins truments onboard t he Herschel Space Observatory 
fHSO. IPilbratt et aLllioTTI ). The dust in those galaxies ab- 
sorbs most UV photons and re-radiates in the far-infrared 
and sub-mm wavelengths. 

In the near-infrared restframe, the K band 2 — 2.3 /xm) 
samples the peak of the emission of the older stars and is 
hence well suited to evaluate th e stellar mass of a galaxy 
(ICole et al.ll200ll : lBen et al.ll2003l ). The weh established link 
between stellar niass an d specific star formation rate (SSFR, 
iRodighiero et al.ll201Ql , for recent evidence with Herschel) 
ensures that the K band is also interesting when investigat- 
ing the SFR. 

One difficulty in combining object catalogues at widely 
different wavelengths lies in deciding which objects are truly 
associated, and which are unrelated foreground/background 
objects. Matching sub-mm sources to objects observed in 
much shorter wavelengths is particularly difficult, because 
the large sub-mm beam sizes, and hence the comparatively 
lower angular resolution, and high confusion noise increase 
the positional uncertainties, which in turn forces us to in- 
crease the radius we employ searching for plausible counter- 
parts. The increased search radius, together with the high 
surface density of objects from optical/near-infrared sur- 
veys, is responsible for the ineffectiveness of a simple closest 
neighbour method. 

A method often applied to previous sub-mm surveys 
consists of first mat ching sub-mm sou r ces to radio or 
mid- infrared sources Jivison et al ] l2QQ7l : iDve et al.1 l2QQ9l : 
iBiggs et al.1 l201lh . before utilising the multi-wavelength 
data that are already available for the radio/mid-infrared 
counterparts, or using the more accurate positions to im- 
prove on the matching technique. This is a dvantageous 
because of the radio /far- i n frared correlation jHelou et al 
19851 : Igainline et al." 2010; "ivison et al." '2010'; Jarvi s et al 
Michalowski, Watson & Hiorth 
the 
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2010; Bour ne et al 
and the 



lower surface density an d the high posi- 
tional accuracy of radio catalogu es (|lvison et al.l l2007l : 
iDve et al.1 [2QQ9I : IPunlop et al.1 12OIOI ). Unfortunately, iden- 
tifying optical/near-infrared counterparts via radio coun- 
terparts is not yet practical for the Hers chel Astrophysica l 
Terahertz Large Area Survey fH-ATLAS. lEales et al.ll2010l ) 
since current radio teles copes cannot del i ver th e required 
area/depth combination: iHardcastle et al.l (|2010l l found 187 
radio sources within the H-ATLAS SDP (Science Demon- 



stration Phase) field, less than 3% of the H-ATLAS sources. 
Radio surveys will dramatically improve in future with SKA 
(Square Kilometre Array) and its precursors LOFAR (LOw 
Frequency ARray, http://www.lofar-uk.org), ASKAP (Aus- 
tralia n Square Kilometr e Array Pathfinder, [ PeBoer et al.l 
20091). and MeerKA T (ISchilizzi. Dewdnev Lazid I2OO8I : 



Johnston et al.|[2008l V b ut not for several years. 



In the mid-infrared, iRoseboom et al.l (|2010h have used 
Spitzer source positions for the source extraction process 
and deblending of HerMES (Herschel Multi- Tiered Extra- 
galactic Survey) sources, taking advantage of the small po- 
sitional uncertainties at 3.6 fim and 24 fim. However, mid- 
infrared data are available for only small patches wi thin the 
H-ATLAS observed fields fsee lBond et alJ (|in prep.!)), using 
WISE (Wide-fie ld Infrared Survey E xplorer, I Wright et al.l 
I2OIOI ) data and iKim et al.1 ([ in DreD.h . using Spitzer-IKAC 
data. 

An alternative approach, adopted here, is to match the 
sub-mm sources directly to a near-infrared catalogue by 
using information on the positional uncertainty probabil- 
ity distribution of the sub-mm sources and the magni- 
tude distribution of the optical/near-infrared objects util- 
ising a likelihood ratio method (LR). This approach uses 
the ratio of the probabilities of a match being the true 
counterpart and being an unrelated background object 
Me Ruite r. Arp W illis "1977"; Prestage & Peacock "l983l: 
Wolst encroft et al.lll986 : Sutherland & S aunders. .1992). We 
describe the likelihood ratio method in detail in section [3Tl 

The paper is structured as follows. Section 2 describes the 
surveys and the data selection. The LR method is discussed 
in detail in section 3, and in section 4 we explain the method 
to obtain our photometric redshifts. Section 5 presents our 
identification results, and in section 6 we compare our results 
with those of the optical matching supplied with the SPIRE 
Phase 1 catalogue. Section 7 summarises our conclusions. 



2 DATA 



2.1 Herschel-ATLAS SPIRE sources 



The Herschel-ATLAS survey (|Eales et al.1 12OIOI ) is a large 
open-time key project carried out with the Herschel Space 
Observatory. The full survey will cover approximately 
550 deg^ of high galactic latitude sky in six patches; the 
survey covers the wavelength range 100 — 500 /xm, provid- 
ing imaging and photometry. Two instruments survey in 5 
passbands, centred on wavelengths 100 and 160 /xm (PACS, 
Poglitsch et alll2010l ) and 250, 350 and 500 /im (SPIRE, 



Griffin et al.ll2010l ). The beams have FWHM (full width at 



half maximum) of respectively 8.7, 13.1, 18.1, 25.2 and 36.9 
arcsec, with 5cr point source sensitivities of 132, 126, 32, 36 
and 45 mJy in the above 5 passb ands. The maps and data 
reduction a re dis cussed in detail in lPascale et alJ (|201lh and 
llbar et al and the so urce catalogue creation is de- 

scribed in lRigbv et aD (l201lh . 

The H-ATLAS fields have been selected for low cirrus 
foreground, and overlap with a number of other existing 
and planned surveys to profit from multi- wavelength data. A 
few important overlapping surve ys are the Sloan Digital Sky 
Survey fSDSS. lYork et al.ll2000l ) in the optical, the Galaxy 
And Mass Assembly (GAMA, [Driver et aD l201lh survey 
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which includes a spectroscopic redsh ift survey of mostly 
SDSS objects (see iBaldrv et alj l20ld for the t arget selec- 
tion) an d the VISTA VIKING imaging survey (jSutherlandl 
lin prep .[) in the near- infrared. 

Here we use the H-ATLAS Phase 1 catalogue of the 
GAMA 9h (G09) field, comprised of 26,269 sources detected 
at 5(7 in the 250//m band. The H-A TLAS catalogue sup - 
plies optical counterparts from SDSS teov os et ahllin prep.! ) 
using a si milar LR technique to that presented here, as dis- 
cussed in lSmith etHI (|2011al ). We find that 22,000 of these 
SPIRE positions are within the region observed by VIKING 
up to late 2010 (an area of approximately 55 deg^, contained 
within RA between 128 and 141 degrees and Dec between 
—2 and +3 degrees), and these comprise our catalogue used 
in the matching hereafter. 



2.2 VISTA VIKING data 

VISTA is a 4m wi de- field telescope at the ESQ Paranal ob- 
servatory in Chile (jEmerson Sutherlandll2010l ). The cam- 
era has 16 near-infrared detectors and an instantaneous 
field of view of 0.6 deg^, and its filter set includes the 
five broad-band filters Z,Y,J,H,Ks with central wavelengths 
0.88 - 2.15 /xm. The VISTA Kilo-degree INfrared Galaxy 
(VIKING, Sutherland et al., in prep.) survey is one of the 
public, large-scale surveys ongoing with VISTA. The sur- 
vey has been planned to cover around 1500 deg^ of extra- 
galactic sky in the above five filters, including one southern 
stripe (including the H-ATLAS SGP stripes), one equato- 
rial strip in the North galactic cap (including the GAMA 
12h and 15h fields) and also the GAMA 9h field. The me- 
dian image quality is ~ 0.9 arcsec, and typical 5cr magni- 
tude limits are J ^ 21.0, Ks ~ 19.2 on the Vega system or 
J ^ 21.9, Ks ^ 21.0 on t he AB system. 

The data processing ([Lewis, Irwin BunclarkI l2010l ) is 
a collaboration of the Cambridge Astronomy Survey Unit 
(CASU) and the Wide Field Astronomy Unit (WFAU) at 
the Royal Observatory in Edinburgh. The data used in this 
paper is from the Vista Science Archive (VSA) produced and 
maintained in Edinburgh, released internally on the 14th 
April 2011; this is our preliminary object catalogue. The 
VSA builds on the WFCAM Science Archive (WSA), pro- 
viding similar access and functionality^ (e.g. image cut-outs, 
SQL queries etc.). Sources are extracted after the merging 
of individual frames and are listed in tables together with 
astrometric and photometric information. 

For our object catalogue, we request a 5cr Kg band detec- 
tion, using (aperture corrected) aperture photometry with a 
diameter of 2^^ We also require a J band detection to exclude 
the large majority of spurious detections (bright star halos, 
satellite trails etc). In addition, we only use sources that are 
primary detections (best source in overlap regions) with er- 
ror bits smaller than 256 (only informational error quality 
conditions), and sources flagged as saturated or noise were 
excluded. This results in 1,376,606 objects in our VIKING 
catalogue from the G09 field. 



^ for a det ailed description of t he functionality and the access 
options, see lHamblv et all (|2008l ) 



2.3 Star-galaxy separation 

The likelihood ratio method we employ to match both cata- 
logues (see section 3) , uses the magnitude distribution of the 
true counterparts which will depend on the morphological 
type of the VIKING objects. The VSA uses a shape param- 
eter calculated from the brightness profile of the objects and 
the point-spread-function (PSF) on each individual detector 
to classify objects as stars and galaxies. The galaxy sample 
is optimised for completeness, leading to a stellar sample 
that is optimised for reliability. We have hence started build- 
ing our stellar sample by using the VSA stellar probability 
> 95%. This puts 464,033 objects firmly into the stellar 
class. 

The remaining objects were then matched to the SDSS 
database, using the nearest neighbour within 2^^, to obtain 
g — i colours for a classification on the J — Kg vs g — i 
colour-colour diagram. This follows the procedure used by 
iBaldrv et al ] (Hoio) to select a highly complete galaxy sam- 
ple for the GAMA input catalogue, in conjunction with a 
SDSS shape parameter. For objects with r < 19.8, we em- 
ploy their prescription for the classification, using a combi- 
nation of shape and colour parameters. However, most of 
ou r objects are much fainter than the r < 19.8 cut used 
bv IBaldrv et al.l (|2010l ), and here the SDSS morphological 
classification is unreliable, as was evident when we used a 
sub-sample with available spectroscopic redshifts. For ob- 
jects with an SDSS counterpart fainter than r > 19.8, we 
then classify via the position on the colour-colour diagram. 
Fig.[T]shows the colours of the VIKING sample, the location 
of our stellar locus and the chosen separation line. 

The remainder of VIKING objects that have not been 
classified above, i.e. objects with VSA stellar probability 
< 0.95 and without an SDSS counterpart, are separated 
as follows: objects with (J — Ks)ab > 0.21 are classified as 
galaxies and those with {J — Ks)ab < —0.34 as stars. The 
logic of this can be seen on the colour-colour diagram: even 
without g — i information, these objects must lie respectively 
above/below the black separation line in Fig.[T] This leaves a 
stripe at intermediate J — Ks colour where the colour classi- 
fication remains ambiguous: just over 12,000 objects fall into 
this category. Here, we look at the VSA shape classification 
again and relax our earlier cut of 95% to 70%, classifying ob- 
jects as stars or galaxies with a cut at pStar = 0.7. Finally, 
we move 573 objects classified as stars to the galaxy class as 
they have confirmed non-stellar redshift of 2; > 0.002 from 
SDSS spectra: those are mostly confirmed QSOs. This clas- 
sification results in a sample of 847,530 galaxies and 529,076 
starfl 

Using this method of classification, QSOs without spectra 
in SDSS are mainly classified as stars, selected by morphol- 



In the catalogue, we include the flag 'sgmode' which indicates 
how we have arrived at the classi fication: 1 VSA star with pStar> 
0.95, 2 uses IBaldrv et al.l (|2010l ) for objects with r < 19.8, with 
slight modifications, 3 colour-colour selection for objects with r > 
19.8, 4 (J — Ks)(^AB) colour selection for objects without SDSS 
counterpart within 2^^, 5 / 6 VSA pStar= 0.7 cut for objects 
without SDSS counterparts and ambiguous in J — Kg colour. The 
sgmode flag is changed by appending a zero to the initial flag 
if the object was moved from the star class to the galaxy class 
having a confirmed non-stellar redshift. 
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Figure 1. — The colour-colour diagram of VIKING objects with 
SDSS counterparts. The red line shows the stellar locus, obtained 
from fitting a quadratic equation in the range [0.3,2.3] to the 
sample with VSA probability of being a star > 0.95. The black 
line is offset +0.3 mag from the locus, representing the adopted 
star/galaxy separation cut. Objects without SDSS identifications 
are classified according to their J — Ks colour only, leaving a 
stripe of —0.34 < (J — Ks)ab < 0.21 where the colour-based 
classification is ambiguous. 



ogy. This is clearly not ideal, because objects in the star 
sample are less likely to be identified as reliable counter- 
parts to the SPIRE sources. A more detailed separation, 
taking QS O properties into acco unt, will be explored in fu- 
ture work (|Hovos et alJIin prep. I ). 



3 APPLYING THE LIKELIHOOD RATIO 
METHOD 

3.1 The Likelihood Ratio method 



One of the earliest approaches jde Ruiter, Arp WillisI 
Il977l : IPrestage k PeacocklfTosi : IWolstencroft et al.lll986l UQ 

the matching of two source catalogues in different wave- 
lengths uses the ratio of two likelihoods: the likelihood that a 
true counterpart is observed at a distance r from the source 
and with magnitude m, and the likelihood that an unrelated 
background object is observed with the same properties: 

P{r,m,id) 
P{r, m, chance) 

The probability P{id\r,m) that an object at distance r and 
with magnitude m is a true counterpart, also called the re- 
liability, is then: 



P{id\r, m) 



L + 1 



(2) 



us ing Bayes' theorem and the th eorem of total probability. 

ISutherland Saunde^ (| 19921 ) extended the likelihood ra- 
tio by incorporating information about other potential coun- 
terparts to one source in the calculation of the reliability, and 
also including additional information (denoted by c) which 
can be, for example, colour information or star/galaxy clas- 
sification. This is especially useful in a situation where the 
matching catalogue has a high surface density so that there 
is a high probability of there being more than one possible 



counterpart. The reliability Rj that the jth candidate for 
one source is the true counterpart is then: 



P{id\r, m, c) 



(3) 



EiLi + il-Q) 

where i runs over all candidates for this source, and Q is 
the probability (for a random source) of finding a genuine 
counterpart above the limiting magnitude of the matching 
survey. In contrast to equation (|2]), this includes informa- 
tion about other candidate counterparts to the source. To 
clarify this distinction, equation [2] is the reliability without 
information about other candidates for the same source (e.g. 
picking one candidate at random from a concatenated list of 
all candidates for a large number of submm sources), while 
equation [3] is the reliability given the set of all candidate 
m atches for one given sou r ce. 

[Sutherland SaundersI (|l992l l define ^(m,c), the proba- 
bility distribution of the true counterparts with magnitude 
m and additional property c, and f{x,y) the probability 
distribution of the source positional errors x, y, normalised 
such that 

/ / f (x^y) dx dy = 1 and / q(m^c)dmdc—Q 

J J —oo J —oo 

(4) 

where Mum is the limiting magnitude in the matching cat- 
alogue. It is usually assumed that the positional errors are 
independent of the magnitude and other additional infor- 
mation, so that P(m, X, 2/, c, id) = P{m,c,id)P{x,y,id). If 
n(m, c) is the surface density of unrelated background ob- 
jects per unit magnitude, the likelihood ratio for any candi- 
date match is then: 

^ ^ P{m,c,x,y,id) q{m,c)f{x,y) 
P{m,c,x,y, chance) n(m^c) 

In practice, the probabilities above have to be estimated 
from the data by fitting simple models. The surface den- 
sity of unrelated background objects is estimated from the 
surface density of objects in the matching catalogue. The 
next two subsections explain how we estimate the distribu- 
tions f{x,y) for the positional errors and q{m) for the true 
counterparts. 

3.2 Positional uncertainties 

We here adopt the simple model that the H- ATLAS source 
positional errors are Gaussian with equal RMS apos in each 
of RA and Dec; then the normalisation condition above re- 
quires 



f{x,y) = f{r) 



exp( 



pos 



(6) 



where r = \J x^ ^ y'^ is the radial positi on difference , and we 
note /(r) has units of (solid angle) ISmith et al.l (|2011al ) 
have estimated the positional errors of > 5cr SPIRE sources, 
using histograms in RA and Dec of the total number of 
SDSS sources within a 50 arcsec box around the SPIRE 
250 /im centres and taken the clustering of the SDSS ob- 
jects into accounl|f|. To be able to use their results, we mea- 
sure the correlation of our VIKING objects to the SDSS 



3 ISmith et aP (|2011al ) fit the sum of the Gaussian positional er- 
rors and the clustering signal, convolved with the Gaussian errors. 
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objects, constructing the corresponding RA and Dec his- 
tograms, with VIKING objects within a box around the 
SDSS positions. The la VIKING position errors are < 0.2'' 
and therefore neghgible compared to the SPIRE errors. We 
hence adopt the weighted mean Icr positional uncert ainty 
of (Tpos = 2.40'' ± 0.09 quoted in ISmith etHI (|2011al ) and 
assume the errors to be symmetric in RA and Dec. 

In theory, the positional uncertainty should depend on 
the SNR ( signal-to-noise ratio ) and the FWHM of the ob- 
servations. Ilvison et al.l (|2007l l deriv the position a l uncer - 
tainty as a = 0.6 x Following ISmith et al.1 (|2011al ). 

we adjust the formula to match our mean positional error 
by inserting a scaling factor of 1.09, so that: 

FWHM 

^vos = 0.655 X (7) 

with the SPIRE mean FWHM= 18.1". We also set a min- 
imum of Gpos = 1" for sources with high SNR, as there 
are limitations to the minimum posi tional accuracy from 
SPIRE and SDSS maps, as discussed in lSmith et aD (l2011al ). 
We adopt a conservative search radius of 10" which would 
include 99.996% of the real counterparts assuming Gaus- 
sian errors; in pract ice, there is evidence for slightly non- 
Gaussian wings (see IHovos et al.l lin prep.! ) , but this radius 
still includes almost all genuine matches. 



3.3 Estimation of q(m^ c) and Q 

We estimate the probability distribution q{m,c) of the true 
counterparts by using the backgro und subtracted sam ple of 
candidate matches, as outlined in ICilied et all (l2003l ). For 
each class (c = galaxies and c = stars), we estimate q{m) 
from the data as follows: 

create a magnitude distribution total{m) of all objects 
within a search radius of 10" around the SPIRE sources 
background subtract total{m) to obtain the so-called 
real{m) distribution 



normalise real{m) so that q{m) 



eal(m) 



real(mi) 



where we sum over bins of magnitude. The normalisation 
factor Qo is an estimate of the probability of finding a coun- 
terpart in the VIKING survey down to the 5cr survey limit, 
Q in equation (|4]13. 

A reasonably accurate value of Qo is important, since this 
enters the reliability formula above, equation [S] Simply es- 
timating Qo via a stacking analysis (summing real{m) and 
dividing by the number of SPIRE sources) is not ideal, since 
source clustering and/or genuine multiple counterparts will 
tend to overestimate Qo by mult i- counting, and therefore 
reliability estimates will be biased high. 

To avoid this multi-counting problem, we decide to es- 
timate 1 — Qo, the fraction of SPIRE sources without 
a VIKING-detected counterpart, hereafter called blanks: 
these will be mostly real sources fainter than the VIKING 



to the resulting histograms. For a more detailed description of the 
derivation of the positional uncertainties, see section 2.1 of their 
paper. 

a derivation of t his formula can be found in the appendix of 
Ilvison et aP (l2007l) 

^ ICiliegi et al .1 (l2003l ) have introduced the constant Qo as the 
value of Q as estimated from the data. 



limit, but also including counterparts outside the search ra- 
dius and spurious SPIRE detections (if any). We start by 
counting the observed blanks to a given search radius r; we 
then need to correct for those sources that have VIKING 
candidate matches, where the match (es) are in fact unre- 
lated to the SPIRE source. The number of true blanks is the 
number of observed blanks, plus the number of true blanks 
that have been matched with a random VIKING object. To 
estimate the latter, we create a catalogue of N (=number 
of SPIRE sources) random positions and cross-match with 
the VIKING catalogue. Hence, defining R as the number of 
blanks at random positions, R the number of random posi- 
tions with a VIKING source within 10" and S the number 
of observed SPIRE blanks, we can calculate the number St 
of true SPIRE blanks as follows: 



1 - R/N R/N 

Dividing by N provides us then with the fraction of the 
SPIRE sources that are true blanks, St/N = S / R, which is 
our estimate for 1 — Qo. Thus, we only need to divide the 
number of SPIRE blanks by the number of random blanks, 
for a given search radius. 

For our default search radius of 10" , we obtain 1 — Qo = 
0.25, or Qo = 0.75. We could use this value in our subsequent 
LR analysis; however, it depends on the value of the search 
radius and the Q, as defined in equation Q, of the VIKING 
catalogue is independent of the search radius. We would like 
to obtain an estimate Qo of Q that is independent of the 
radius and so repeat the above procedure for radii in the 
range 1-15". The values we obtain for the fraction of true 
blanks are shown in Fig. [2] as black points. 

We then model the dependence of the true blanks on the 
search radius as follows: a SPIRE blank at radius r is a 
source whose counterpart is either fainter than the VIKING 
limit, or lies outside the search radius, or both. The for- 
mer probability is 1 — Qo, this is the first term in equa- 
tion . The probability of the counterpart to reside outside 
the search radius can be calculated from the positional error 
distribution /(r), leading to our second term in equation Q. 
The third term follows if we assume that both possibilities 
are independent of each other, and using the standard prob- 
ability result P{A or B) = P{A) + P{B) - P{A and B) . Our 
model for the dependence of the true blanks on the search 
radius r is then: 

(1 - Qo) + (1 - F{r)) - (1 - Qo)(l - F{r)) = 1 - QoF(r) 

where F{r) — I P{r')dr' = 1 — e 2^ 
Jo 

and P(r) = 2nrf{r) 

(9) 

Fitting this model to the data, we obtain Qo = 0.73 ± 0.03 
as our best-fitting value. Fig. [2] shows the best fit model as 
the black line and the black filled circles as our data points 
{S/R for each radius). 

The model underestimates the number of SPIRE blanks 
in the data in the range A" < r < ^" and overestimates 
the data for r > 10" . This might show some evidence for 
clustering of the VIKING objects which we have not explic- 
itly considered, but which is accounted for in the value of 
the mean positional uncertainty a^os = 2.4^^ bv lSmith et al.l 
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Figure 2. — Illustration of the procedure to estimate 1 — Qo- 
The red crosses show the fraction of blank SPIRE sources within 
search radius r, while the blue stars show the same fraction 
for random positions. The black filled circles represent the data 
points obtained from dividing the number of blank SPIRE posi- 
tions by the number of blank random positions, our estimate for 
(1 — QQ)F{r). The black line represents the best fit to the model 
(1 - Qo)F{r), see equation g)), with Q = 0.7342 ± 0.0257. 



(l2011al l: they convolve their model with the clustering sig- 
nal of the SDSS sources to obtain the value we adopt here. 
It might also demonstrate that the Gaussian approxima- 
tion for the SPIRE positional errors, our equation dH}, is 
not entirely accurate. This is also evident when investigat- 
ing histograms of distances of VIKING and SDSS objects to 
SPIRE sources. We see slightly higher numbers of objects 
at distances of around 10" than expected if we assume a 
Gaussian error distribu tion. This assumption is examined 
and will be discussed in lHovos et al.l (|in prep.[ ). 

Our fitted value of Qo = 0.73 is consistent with the value 
of Qo = 0.75 from the datapoint at our search radius of 10^^ 
and is more conservative. We hence adopt the fitted value 
Qo = 0.73 for our subsequent LR analysis. 

A drawback of this approach is that we cannot separate 
our data into stars and galaxies, as we use blank fields. We 
estimate the val ue of Qo for s t ars, de noted Qos in applying 
equation (10) of ISmith et al ] (l2011al ). using a background 
subtracted sample of possible matches within 10", yielding 
a value of 0.004. This is very small indeed and shows how 
unlikely it is that stars are detected with SPIRE. For our 
LR analysis we then adopt Qo^ =0.72 and Qos = 0.01, the 
values for Qo used in the normalisation of q{m) for galaxies 
and stars respectively. 



3.4 Probability of mis-identifying a true 
counterpart 

Given the above model, we now estimate 
P(wrongID\mtrue), defined to be the probability that 
a true counterpart with a given VIKING Ks magnitude 
TTitrue is uot the bcst candidate using our LR method. This 
situation occurs, if the true counterpart has a likelihood 
ratio value L = Li and there exists a chance match with 
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Figure 3. — The solid red line shows the distribution D(L) of the 
highest L values of candidate matches to > 10^ random positions 
in the G09 field. The peaks at lower and higher L values are 
due to star and galaxy matches respectively. The black histogram 
represents the LR values of all candidate matches to the SPIRE 
positions. The filled red histogram represents the reliable matches 
only. 

L > Li for the same source. Hence: 

rLmax 

P{wrongID\mtrue) ^ / P{Li\mtrue) 

Jo (10) 

X P(chance > Li) dLi 

The probability P (chance > Li) that there exists a 
chance match with LR > Li for one SPIRE source can be 
estimated through simulations. The steps of the procedure 
we have used are as follows: 

• Create N random positions on an area common to both 
VIKING and H- ATLAS in the G09 field 

• For each random position, calculate the likelihood ratio L 
for random matches, if any, in the VIKING catalogue. 

• Create the distribution D(L) of the highest likelihood ratio 
value for each random position. 

From the probability distribution D(L) of the highest LR 
values for candidate matches to random positions, we can 
calculate the probability that a given source has a chance 
candidate above Li by chance: 

1 pLmax 

P{chance > Li) = — D{L) dL (11) 

J Li 

A similar method has been employed bv lDve et al.l (|2009l ) 
to calculate the probability of radio associations to the 
BLAST (Balloon-borne Large Aperture Submillimetre Tele- 
scope) submm sources. Fig. [3] shows the resulting distri- 
bution D{L) from > 10^ random positions, together with 
the histogram of the likelihood ratio values of the candidate 
matches to the real SPIRE positions and of the subset of re- 
liable counterparts (e.g. objects with a probability of > 80% 
of being the true counterpart, see section [5T]) . 

The probability P{Li\mtrue), the first factor in the in- 
tegral in equation (p^ . that a true counterpart with Ks 
magnitude mtrue acquires the likelihood ratio Li can be 
calculated analytically from the probability distribution 
Pr{r) = 27vrf{r) of the offsets. For a given m = mtrue 



Herschel- ATLAS: VIKING counterparts 7 



we have fi = f(r) = Li x n{mtrue) / q{mtrue) and hence 
Pfifi)- We can find the probability distribu- 
tion Pf{f) from Pr{r) by performing a variable transforma- 
tion: 



Pfif) = -Pr{r{f)) X 



dr 



2 A 



(12) 



This is a surprisingly simple result, i.e. for Gaussian er- 
rors, Pfif) has a uniform distribution between and its 
maximum value l/27r cr^. We would like to note though that 
the positional uncertainty depends on the value of the SNR. 
For simplicity, we have used the value of apos = 2A" ^ al- 
lowing for most of the SPIRE sources to have a SNR value 
close to 5. A more detailed analysis could take into account 
the probability distribution of the SNR values of all SPIRE 
sources. 

Having calculated the probability of a true counterpart 
with assumed magnitude not being the best candi- 

date in our LR method, we can then calculate P{wrongID), 
the same probability integrated over the model magnitude 
distribution of the true counterparts: 



P{wrongID) 



jMurn p(^yj'pQjigJD^fji^ q(m) dm 
S^'^-^ q{m) dm 



(13) 



Using the distribution D[L) from the random positions, our 
VIKING sample as described in section 12.21 to derive the 
magnitude distributions q{m) and n{m), as well as (T-pos — 
2A" ^ we obtain P(wrongID) — 0.0493. We hence expect to 
mis-identify around ~ 5% of the true VIKING counterparts 
to the SPIRE sources. 

This is likely to overestimate the true value because we 
have not taken the individual SNR values of the SPIRE 
sources into account. This performance measure of the LR 
method is compared in section 15.11 with the false identifica- 
tion rate for reliable counterparts. 



4 PHOTOMETRIC REDSHIFTS 

The submm wavebands benefit from negative k-correction 
with the effect that objects can b e detected at least out 
to z > 5 in the 870 /xm band (jBlain Longa 3 Il993l : 
IChapin et ahlboill ) . Selecting sources at 250 fim reduces the 
effectiveness of the negative k-correction, but models still 
predict a sign ificant fraction of the SPIRE sources t o reside 
at 2; > 1 fe.g. lAmblard et al.ll2010l : iLapi et al.ll201lh . 

A significant fraction of our VIKING counterparts 
28%) have spectroscopic redshift from the GAMA and SDSS 
redshift surveys. For the remainder, we will obtain photo- 
metric redshifts by combining the near-infrared photometry 
with optical photometric information from the SDSS sur- 
vey. Objects with spectroscopic redshift from the GAMA 
survey are then employed as a training set to estimate the 
photometric redshifts using neural networks, see below. 

SDSS matches are found for - 68% of the VIKING can- 
didates within 2" by performing a simple nearest neighbour 
match. The remainder of the VIKING objects are too faint 
to be detected in the SDSS survey. For near-complete visi- 
ble detections we have to wait for the VST KIDS survejQ, 
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which will observe the VIKING areas down to 24.8, 25.4, 
25.2 and 24.2 (lOcr, AB) in u\g\r\i bands respectively, 
thus giving detections in at least r\i bands for nearly all 
VIKING objects. 

With a search radius of 2" we are unlikely to identify 
a wrong SDSS counterpart. Using the surface density of 
the SDSS catalogu e employed for the optical matching by 
ISmith et al.l (|2011al ). we estimate around 2% (170 objects) of 
our subsample of reliable VIKING matches to have a wrong 
SDSS identification within 2^^ 

For photometric redshifts, we use a neural network 
method in which the photometry of objects with available 
spectroscopic redshifts provides the sample to train a net- 
work. Once trained, the network is then used to obtain pho- 
tometric redshifts from the photometry of objects without 
reliable redshift information. The photometry in the differ- 
ent bands employed can differ, for instance we are using 
SDSS modelmags and VIKING Vega aperture magnitudes. 
In contrast, a template fitting method where the photom- 
etry is compared to the expected photometry from a class 
of different (empirical and/or theoretical) SEDs (Spectral 
Energy Distribution), needs to use very carefully calibrated 
photometry in all bands. The s oftware used to obtain the 
photometric redshifts is ANNz (jCollister Lahavl [20041 ) , a 
publicly available product. We train a committee of 3 net- 
works for each possible photometric band combination, rec- 
ommended to minimise the network variance. The output of 
ANNz is the photometric redshift for each object together 
with a redshift error estimate which takes the photometric 
errors in each band into account. 

We use 30,697 objects with photometry from SDSS 
and VIKING and spectroscopic redshifts from GAMA in 
our training set. The median spectroscopic redshift of the 
GAMA sample IS Zspec 

= 0.211. The median redshift of 
the VIKING counterparts is expected to be higher and 
the training set should reflect this. Currently, there are no 
deeper spectroscopic surveys available that overlap with the 
VIKING area. There are though some deeper spectroscopic 
redshifts fro m surveys t hat a r e within the UKID SS LAS 
area (.Lawrence et al.1 12OIOI : iHewitt et all I2OO6I ). There- 
fore, we undertake a comparison of WFCAM and VISTA 
photometry, to be able to u se deeper photom etric informa- 
tion from the zCO SMOd!! (iLillv et ahllioOTl ) and DEEPfl 
(iDavis et al.ll2003l ) surveys. We match VIKING objects in 
a 4 deg^ area in the H- ATLAS SDP field to objects in the 
UKIDSS LAS and obtain the mean and standard error for 
the difference in magnitudes (apermag3) in each of the bands 
YJHK. We are then able to use the LAS photometry for 
these samples by subtracting this mean from the LAS mag- 
nitudes and adding the standard error of the difference in 
quadrature to the photometric error. This allows us to use 
the deeper photometry and spectroscopic redshifts of zCos- 
mos (1530 objects) and DEEP2 (238 objects), which overlap 



^ UKIDSS LAS is carried out with the UKIRT WFCAM, and 
images an area of ~ 2500 deg2 in the YJHK filters to a depth 
K=18A 

s zCOSMOS is a redshift survey carried out on the ESQ VLT 
with the VIMOS spectrograph on the COSMOS field 
^ DEEP2 is a redshift survey carried out with the Keck telescopes 
with a pre-selected redshift range of 0.75 — 1.4. 
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with the UKIDSS LAS survey (but not with VIKING), as a 
training set. 

Adding all our subsamples together, we obtain an over- 
all training catalogue with 32,465 spectroscopic redshifts. 
We then compile a photometric catalogue for all VIKING 
candidate matches with ugriz photometry (modelmags) and 
YJHK photometry (aperture magnitudes), where present. 
We attempt photometric redshifts where we have at least 2 
infrared bands with good photometry. 

Using the deeper spectroscopic redshifts from zCOSMOS 
and DEEP2 forces us to exclude the VIKING Z band in 
the training catalogue. Investigating how this influences our 
photometric redshifts, we create a second training catalogue, 
this time using just the GAMA subset (95% of the first 
training catalogue), so that we can include the VIKING Z 
band. The scatter in the differences between photometric 
and spectroscopic redshifts for objects with a spectroscopic 
redshift is slightly lower for this second training set, but this 
is to be expected because we compare only the lower redshift 
end. Assuming that the inclusion of deeper redshifts into the 
training set reflects the true VIKING redshift distribution 
better, we adopt the photometric redshifts from our first 
training catalogue, the GAMA-zC0SM0S-DEEP2 training 
set. 



5 RESULTS 

5.1 Ks band matching 

There are 22,000 SPIRE sources within the sky area corre- 
sponding to the G09 VIKING object catalogue. Of those, 
18,989 sources have at least one possible match within a 
search radius of 10^^ in the VIKING Ks band selected cat- 
alogue, with a total of 35,800 candidate matches; of which 
30,659 are classified as galaxies (85.6%) and 5,141 are clas- 
sified as stars (14.4%), as described in section [231 Table [T] 
shows the number of SPIRE sources matched as a function 
of the number of candidate matches found per position. 

There are 11,294 SPIRE sources with a best VIKING 
counterpart with a reliability of > 80% (11,282 galaxies/ 12 
stars). This means we were able to match ^ 51% of the 
SPIRE sources with a high reliability. We will refer to the set 
of matches with R > 0.8 as "reliable" hereafter (as we show 
later, the mean reliability of this set is near 0.95). Fig. [H 
shows the magnitude dependent distribution q{m)/n{m) for 
galaxies and stars used in the calculation of the LR values. 
Fig. [5] and [6] show the magnitude distributions involved, of 
the background n(m), of the possible matches total (m) and 
of the background subtracted sample real{m), in estimating 
q{m). Also shown in the figures is the magnitude distribution 
for reliable counterparts. The reliable matches show a lower 
fraction of fainter counterparts compared to all candidate 
matches, representing the steep increase of fainter objects 
in the background number counts, causing lower reliability 
values for fainter matches. 

We can estimate a false ID rate by summing up the com- 
plement of reliability values of the reliable matches: 

NfalselD = ^(1 - R^) = 469.25 (14) 

i 

corresponding to a mean reliability of 0.958 and a false ID 
rate of 4.2%. 



100000 




100 : 

10 r 



1 I — ^ ^ ^ ^ ^ ^ ^ ^ 

12 13 14 15 16 17 18 19 

Ks (VEGA) 

Figure 4. — The magnitude dependent part q(m)/n{m) of the 
LR analysis, calculated from the data as described in section [3?3] 
The black line indicates the galaxy candidates, the red line shows 
the stellar candidates. The values for the bright bins for the galaxy 
distribution were extrapolated from the first bin that included 
more than 10 galaxies, at m = 13.4. 
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Figure 5. — The magnitude distributions involved in estimat- 
ing q(m) for galaxies. Here n(m) (solid black) is the distribution 
of the background objects, calculated from the whole VIKING 
G09 catalogue, as described in section [231 and total{m) (dashed 
blue) is the magnitude distribution of all possible matches within 
10" . real(m) (solid red) is the background subtracted distribu- 
tion as described in section I3.3l and is significantly brighter than 
the background. The light grey shaded histogram represents the 
magnitude distribution of the reliable matches. 

In the appendix, we show cut-outs around SPIRE sources 
and corresponding reliable VIKING sources for 16 positions, 
drawn at random from our reliable sample, in magnitude 
bins of 0.4 in the range of Ks = 13.4 - 19.8. 

5.2 VIKING and SPIRE colours 

Figure [71 shows the ZJK colour-colour diagram of the 
10,121 reliable counterparts (red) with 5cr detection in all 
3 VIKING bands. Golours from randomly selected back- 
ground galaxies are depicted in grey. The redshift evolution 




Figure 6. — The magnitude distributions involved in estimating 
q{m) for stars. The colour scheme is as in Fig. [S] 



Table 1. The distribution of the number of VIKING Ks band 
sources within 10 arcsec of the 250 /im SPIRE positions. Of the 
22,000 sources in the VIKING area, 8118 have only one possible 
match within 10 arcsec and 59.7% of these are determined to 
be reliable. This emphasises the difference from a simple nearest 
neighbour match. 



N (matches) 


N(SPIRE) 


N (reliable) 


% 





3011 
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8118 


4851 


59.76 
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6619 


4040 


61.04 
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2968 


1710 


57.61 
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968 


529 


54.65 


5 


241 


128 


53.11 


6 


63 


31 


49.21 
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11 


5 


45.45 
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1 





00.00 


Totals 


22,000 


11,294 





of the s ub-mm selected mean galaxy template o f ISmith et alJ 
l!2011bh is shown in green. The template has been artificially 
redshifted between 2; = — 1.5 in intervals of dz = 0.1 and 
colours have been computed by integrating the product of 
the template SED with the VISTA response functions at 
each redshift interval. A small deviation from the Vega sys- 
tem is present in the V ISTA Z-band and a measured offset 
([Pindlav et al.|[in pressh has been added to the colours com- 
puted here to reflect this. 

The median Z — J and J — Ks colours for the reliable 
matches are 0.97 and 1.64; for the background objects the 
median colours are 0.89 and 1.54. Performing a two-sided 
K-S test on the J — Ks and the Z — J colours for reliable 
matches and the background objects enables us to reject at 
a significance of 99.96% that the two populations are drawn 
from the same distribution for either colour. Thus, we find 
evidence that the reliable matches to the SPIRE sources are 
slightly, but significantly redder than the population of all 
VIKING objects in the G09 field. 

Within the search radius of 10^^, 3011 sources have no 5cr 
Ks VIKING candidate, i.e. they are almost certainly fainter 
than the VIKING limit. Fig. [8] shows 5'25o/*S'350 histograms 



Figure 7. — VIKING colour-colour diagram of reliable coun- 
terparts (red) and background galaxies (grey) with 5cr detections 
in Z,J and Ks. The horizontal and vertical histograms show the 
distributions in J — Ks and Z — J respectively. The green points 
connected by a green line represent the evolutionary track of a 
typical H-ATLAS galaxy SED for z < 0.35 from Smith et al. 
(2011b) and calculated by Findlay et al. (in press). 

separately for the sources with a reliable match (best R > 
0.8), sources with non-reliable match(es) (best R < 0.8), 
and for sources that are blank in VIKING (black, red and 
shaded red respectively). The blank sources (median=1.01) 
have distinctively redder colours than the sources with reli- 
able matches (median=1.32), suggesting that they reside at 
higher redshifts with the peak of the dust emission moving 
to longer wavelengths. The colours of sources with unreli- 
able matches (median=l.ll) lie in between the other two 
populations, indicating that they might be composed from 
members of both populations. From our value of Qo, we ex- 
pect around 60% of the unreliably matched sources to have 
a true counterpart, but for which we do not have a high 
enough reliability, with the remai ning 40% be i n^ ma tched 
to unrelated background objects. ISmith et al.l (|2011al l find 
a similar trend in the distribution of the 5'25o/*S'350 colour 
of SPIRE sources in the SDP field matched to the SDSS 
r-band catalogue for the 3 different populations. 

We also show a SPIRE colour-colour diagram of the 
sources in Fig. [9] The colours of the reliably matched and 
the bla nk sources are very similar to those in ISmith et alJ 
(I2OI lap, their fig. 9. This figure can also be compared to fig. 
1 in Amblard et al.l (|2010l V We add the evoluti onary tracks 
of two t emplates: a low-z template compiled bv lSmith et al.l 
(|2011bh from optical counterparts to SPIRE sources out 
to z ^ 0. 3 5 (bl ue line) and the sub-mm template from 
iLapi et al.1 (l201lh ( green line), thought to be appropriate 
for high-z H-ATLAS sources at z > 1.2. From the tracks of 
the two SED templates we can again suggest that the blank 
sources lie at higher redshifts in general than the sources 
with reliable matches. Assuming that the H-ATLAS sources 
are comprised of two distinct populations, see section 15.51 
a lower redshift population with mainly normal galaxies 
(with a much higher star- format ion rate than local normal 
galaxies), and a higher redshift population of dusty sub-mm 
galaxies (likely to be giant proto-spheroidal galaxies in the 
process of forming most of their stars), our blank sources 
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Figure 8. — Distribution of the SPIRE S'25o/'S'35o colour. Black: 
sources with reliable matches. Red: sources that have unreliable 
matches only. Red shaded: blank sources. The blank sources show 
distinctively redder colours than the other two populations, sug- 
gesting they are at higher redshifts. 
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Figure 9. — Colour-colour diagram of SPIRE objects with 5cr 
detection in 250/im and 350/im and 3cr detection in 500/im. The 
dark grey points represent SPIRE sources that have reliable coun- 
terparts, the red points show the colours of blank SPIRE sources. 
The blue points connected by the blue line represen t the colours 
of a typ ical H-ATLAS galaxy SED for z < 0.35 from lSmith et al.l 
(|2011bl '). The dots are at an interval of z = 0.5, starting at 
2; = on the right hand side. The green points conne cted by the 
green line represent the colours of the sub-mm SED in lLapi et al.l 
(|2Qllh . The points are at an interval of z = 1, starting with z = 
on the right. 

could represent a mixture of the former at z > 0.7 and the 
latter at 2; > 1.5. 



5.3 Multiple counterparts 

The LR method assumes that there is only one true counter- 
part to each source, and assigns reliabilities self-consistently 
based on this, so that the sum of reliabilities cannot exceed 1. 
Thus, if more than one counterpart with > 0.2 is present, 
we will not find a reliable match. If individual reliabilities 



add up to our threshold of 0.8, we could assume that these 
candidate matches are all associated with the sources, either 
through confusion or in a real physical sense (i.e. merging 
galaxies) . 

In our results, we find 1444 SPIRE sources that fulfill 
above the criteria and potentially have multiple true coun- 
terparts in VIKING. Most of those will have additional close 
chance objects within the search radius, denying a reliable 
identification, but some might be genuine mergers or consti- 
tute members of the same cluster. We can rule out a chance 
match by comparing the redshifts of all possible matches to 
one SPIRE source. Checking for available redshifts in the 
GAMA and SDSS spectroscopic redshift databases, we find 
matches to 37 sources whose redshifts are within 5% of each 
other. The mean redshift difference is 0.0011 with a maxi- 
mum difference of Az — 0.0187. Those could be either merg- 
ing galaxies or members of the same cluster. 

We use our photometric redshifts, see section 15.51 to se- 
lect further candidates. We account for the higher errors in 
the photometric redshifts by allowing a redshift difference 
of 10% and also compare photometric redshifts of possible 
matches for which we do not have a spectroscopic redshift 
with spectroscopic redshifts of other candidates to the same 
source. We find 602 further sources where the SPIRE flux 
potentially originates from an interacting system or from 
galaxies within the same cluster. Fig. [10] shows a VIKING 
Ks image cutout around one of the those sources, H ATLAS 
J091017. 1-005538. Due to the uncertainties in the photo- 
metric redshifts, those sources can only be regarded as can- 
didates and need further investigation to be confirmed as 
interacting systems or as members of the same cluster. 

It would be interesting to confirm how many sources def- 
initely do not have physically related multiple counterparts 
and are just unreliable matches, but this is difficult due to 
the sparsity of available spectroscopic redshifts and the un- 
certainties on the photometric redshifts. However, we can es- 
timate the number of reliable identifications we are missing 
due to potentially multiple counterparts. From Table [1] we 
can see that the identification rate for reliable counterparts 
is approximately 60% without the presence of additional po- 
tential matches. We do not see a decrease in the identifica- 
tion rate for sources with two possible matches, suggesting 
that the true number of merging galaxy pairs is indeed low. 
For sources with higher numbers of possible matches, we 
have an increased possibility of having observed a galaxy 
cluster and so the identification rate for reliable matches 
falls. For instance, sources with 4 possible matches have an 
identification rate of around 55%, suggesting that we miss 
5% of the reliable matches, equivalent to 50 sources. Adding 
up the missed reliable matches of all sources with more than 
2 possible counterparts, suggests that we are missing around 
150 reliable VIKING counterparts due to additional matches 
within our search radius. This is a small number indeed, 
only around 10% of the number on our list of candidates 
for true multiple counterparts. There is evidence, from ob- 
servation and simulation, for the merger rate to evolve with 
reds hift and to p eak z ^ 1.2 (iBell et al.ll2006l : lLotdl2"QQ7l : 
Rvan et al.ll2008l l , beyond the redshift of most of our reliable 
counterparts. Hence, we might miss a substantial fraction of 
mergers not because we find multiple candidates but rather 
because they hide in the fraction of blank SPIRE sources. 
This implies that our candidate list comprises mostly either 
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Figure 10. — VIKING Ks image of HATLAS J091017.1-005538, 
15" on the side. The black cross indicates the SPIRE position. The 
blue circled object has reliability R = 0.74, the red circled object 
has R = 0.25. The third candidate has negligible reliability, but 
could still be part of the interacting system. 



chance alignments of galaxies or cores of clusters of galaxies 
resulting in confusion when observed with SPIRE. 



5.4 Stellar matches 

The far- infrared /sub- mm mainly detects cold, dusty objects. 
It is unlikely that stars are detected with SPIRE, unless they 
are post-AGB, s hrouded in dust or have debris disks (e.g. 
[Thompson et al.l^2010l We have matched 12 SPIRE sources 
reliably to point-like objects. Their location on the J — Ks 
vs g — i colour- colour diagram is displayed in Fig. [TT1 

HATLAS J090450.4-014525 (884 in Fig. [HI) displays 
galaxy- like col ours and is listed as a Q SO in the quasar 
catalogue of iVeron-Cettv VeronI (|2010l ). All other objects 
are consistent with having star-like c olours. 

Tw o blazars were identified by I Gonzales- Nuevo et al.l 
(|2010l ) in the SDP field from cross-matching to radio ob- 
servations. We match H- ATLAS J090910. 1+012135 reliably 
to a point-like VIKING object (sgmode= 10, i.e. point-like 
object classified as a galaxy on the basis of a non- stellar 
spectroscopic redshift). H- ATLAS J090940.3+020000 is not 
matched reliably within our search radius of 10". We find 
one possible VIKING counterpart within 10" of this SPIRE 
source: a bright point-like object {Ks = 15.03, sgmode= 10), 
lying at a distance of nearly 9" from the SPIRE position and 
therefore obtaining a low reliability. However, it lies within 
1" of the known blazar PKS 0907 +022, the object identified 
as the blazar counterpart to H - ATLAS J090940.3+020000 
bv lGonzales-Nuevo et aD (|2010l ). Despite the low reliability, 
there is hence evidence that our VIKING object is the coun- 
terpart to H- ATLAS J090940.3+020000. The colours of both 
VIKING objects are shown in Fig. [TT]as blue dots. Blazars 
in the H-ATLAS Phase 1 fields are currently investigated 
by a team led by Marcos Lopez-Caniego with 14 candidates 
identified so far. 

For the brighter stars, it is a possibility that the measured 




Figure 11. — Reliable stellar matches, represented by the red 
dots. For clarity we have labeled the dots with row numbers from 
the SPIRE catalogue and not with the full HATLAS names. HAT- 
LAS J090450. 5-014525 (884) has galaxy-like colours and was orig- 
inally selected to be a QSO target of the Sloan spectroscopic sur- 
vey but has subsequently been rejected as a target. It is listed 
as a QSO in the Quasar and Active Galactic Nuclei catalogue by 
IVeron-Cettv vironi (|2Q1G| ) with redshift z = 1.005. The colours 
of the other objects are consistent with being star-like. The ex- 
treme J-Ks = -2.67 value of HATLAS J091233. 9-004549 is not 
shown in the diagram and could be due to saturation. The blue 
dots represent the colours of the VIKING counterparts to the two 
blazars found in the SDP field, 1 =H- ATLAS J090910. 1+012135, 
2 =H-ATLAS J090940. 3+020000 



SPIRE flux originates from a galaxy that is too faint in the 
near-infrared to be detected by VIKING, or obliterated by 
the star, and the star is a chance projection. The reliabilities 
of reliable star matches are on average lower (R = 0.85) than 
for the reliable galaxy matches, due to the lower Qo value 
and the lower values of q{m) (see Fig. |4]). 



5.5 Photometric redshift distribution 

We have measured the photometric redshifts for all of our 
11,294 reliable VIKING counterparts, 8750 of which have 
SDSS matches within 2", as described in section 4. Spec- 
troscopic redshifts exist for 3147 of the reliable VIKING 
counterparts which allows us to evaluate the accuracy of 
the photometric redshifts, see Fig. [T2l The spectroscopic 
redshifts are taken from either the GAMA or the SDSS red- 
shift survey. Where a redshift exists in both surveys for a 
VIKING object, we have used the GAMA redshift (quality 
flag Q > 3 only). Excluding 31 confirmed QSO with z > 1, 
the scatter of the difference between our photometric red- 
shifts and the spectroscopic redshifts is a = 0.0353. This 
reduces to a = 0.0259 for the normalised redshift distribu- 
tion Az/{1 + Zspec)- 

Also visible from Fig. [12] is the tendency for the photo- 
metric redshift to underestimate the redshift at higher val- 
ues of z] for z > 0.3, the bias amounts to —0.0889. The 
systematic underestimation is also found when comparing 
spectroscopic redshifts with the photometric redshifts from 
the H-ATLAS Phase 1 catalogue which used optical and 
near- infrared photo metry from SDSS and UKIDSS-LAS, see 
ISmith et al1(|2011al ). They have employed a similar training 
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Figure 12. — Upper panel: photometric vs. spectroscopic redshift 
for 3147 reliable VIKING counterparts with spectroscopic redshift 
from either GAM A or SDSS. Lower panel: scatter of the difference 
\zspec — Zphotl for the same 3147 objects in redshift bins of Sz = 
0.05. 



set which suggest that we face an issue with the representa- 
tiveness of our training set. The reason for the bias seems less 
likely to be a lack of spectroscopic redshift s at z > 0.3 but 
could rather be related to a difference in the colour distribu- 
tion of galaxies in the training set and our VIKING galaxies. 
Clearly, more work is needed to investigate the reasons for 
the bias at z > 0.3 and, crucially, to assemble a more rep- 
resentative training set which is outside the scope of this 
paper. 

The redshift differences Zphot — Zspec of the 31 confirmed 
QSO with z > 1 are considerably worse, as the training set 
includes few high-z QSO and also the near power-law spec- 
tra of QSOs means that QSO photo-z estimates are much 
worse than for galaxies. At the time, we cannot compare 
photometric redshifts for VIKING objects without a SDSS 
counterpart, calculated from just the near-infrared photom- 
etry, due to the lack of spectroscopic redshifts. In general, 
they have higher photometric redshifts, as can be seen from 

Fig.ini 

Certainly, we are missing more reliable identifications at 
the higher redshift end than at lower redshifts. The lower 
panel in Fig. 1 141 shows a comparison of the expected redshift 
distribution, real{zphot) of the VIKING counterparts to the 
redshift distribution of our reliable counterparts {R > 0.8). 
The expected distribution real{zphot) can be calculated in 
a similar way as the magnitude distributions: from the 
total {z phot) ^ the photometric redshift distribution of the 
VIKING objects within our search radius of 10^^ and the 
background distribution n{zphot), both shown in the upper 
panel of Fig. 1141 Table [2] shows the fraction of reliable to ex- 
pected counterparts per redshift bin, i.e. the completeness 
of our photomet r ic red shift sample. 

iDunloD et all tOld ) have calculated very accurate pho- 
tometric redshifts for counterparts to 20 bright BLAST 
sources, using a wide range of deep mult i- wavelength data in 
the GOODS-North field. BLAST used a proto-type SPIRE 
camera and the BLAST sources were selected down to 
36mJy in 250 /xm. We compare our reliable sample to their 
data on a Kg vs Zphot in Fig. [151 This plot suggests that 



Figure 13. — Redshift distributions for the reliable VIKING 
matches. The filled boxes represent the spectroscopic redshifts for 
3147 of the 11,294 rehable VIKING matches with either GAMA or 
SDSS spectroscopic redshift with a median redshift of z = 0.199. 
The excess peak of all photo-z at z ~ 0.5 (black histogram) is 
formed of redshifts to VIKING objects without optical photome- 
try. 
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Figure 14. — Upper panel: the dashed red histogram shows the 
photometric redshift distribution of VIKING background objects, 
the black line represents the photometric redshift distribution of 
all VIKING objects within our search radius (total(zphot))- Lower 
panel: the expected photometric redshift distribution of the true 
VIKING counterparts (real(zphot^ black line), calculated by sub- 
tracting the background from total(zphot- The filled histogram 
represents spectroscopic reshifts and the red dashed histogram is 
the distribution for our reliable counterparts. 



we are fairly complete out to redshift z ^ 1 and are miss- 
ing higher redshift counterparts due to the VIKING survey 
limit. 

Work carried out bv lAmblard et al ] (|2010h using SPIRE 
and PACS colours, suggests that the SPIRE source red- 
shift distribution might be bimodal, formed by a low red- 
shift population of spirals and a high redshift population o f 
starburst galaxies peaking at 2; 2. IPariush et al.1 (l201lh 
found that most of the H- ATLAS low-z galaxies are com- 
prised of blue/star- forming galaxies with some highly dusty. 
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Figure 15. — Ks magnitudes (Vega) versus photometric red- 
shift. Grey dots: our rehable VIKING counterparts. Red points: 
BLAST sources of lDunloD et~all tOldi ) . This suggests that we are 
fairly complete out to 2; ~ 1. 



Table 2. Photometric redshift completeness, calculated from 
a background subtracted sample of possible counterparts, 
real(zphot)- The errors assume that Poisson errors in the real 
distribution are dominant. 



^phot 


Completeness(%) 


O'comp 


0.0-0.1 


117.8 


5.8 


0.1-0.2 


108.6 


3.8 


0.2-0.3 


87 


2.9 


0.3-0.4 


65.9 


1.4 


0.4-0.5 


50.6 


1.3 


0.5-0.6 


46.3 


0.9 


0.6-0.7 


59.3 


1.8 


0.7-0.8 


56.1 


2.7 


0.8-0.9 


55.8 


3.6 


0.9-1.0 


61 


5.6 


1.0-1.1 


55 


5.6 



red spirals. Evidence for this bi modality in the redsh ift dis- 
tribution can also be found in iMaddox et al usms; 
the angular correlation fu nction of 250, 3 50 an d 500 /xm se- 
lected SPIRE sources, in iDunloD et atl (|2010l ). albeit with 
a sample of just 20 BLAST 2b0/j.m s ources, and from the - 
oretical models (iLagache et al. I I2OO3I : iNegrello et ahllioOTl ). 
ISmith et all (|2011al ) discuss this in more detail in their sec- 
tion 3.3 on the photometric redshift distribution. 

In this work, we find no evidence for a substantial popu- 
lation at higher redshifts than z ^ 1 which is explained by 
the photometric depths of our VIKING catalogue to about 
z < 1. From our value of Qo, we expect about 27% of the 
SPIRE sources to be too faint to be detected in VIKING and 
hence very likely to be at higher redshifts, see section 15.21 
Taking into account the possible underestimation of our pho- 
tometric redshifts ai z > 0.3, we find that ^ 2% of the 
H- ATLAS 250 /xm sources with a reliable counterpart lie 
at z > 1. At least a similar fraction of sources without a 
reliable match and expected to have a true counterpart in 
VIKING, should obtain ^ > 1. We hence expect > 30% of 
our H- ATLAS sources to be found at z > 1. We compare this 



with the redshift distributions found for BLAST sources and 
fo r sources detected by SPI RE in the GOODS-N orth field. 

iDunloD et all (|20ld ) and lChapin et al.1 (l201lh . using 250 
fim BLAST sources, both find that ~ 50% of their sources 
lie at 2: > 1 from a variety of photometric redshifts, 
even though the sh apes of the distributions differ, with 
iDunlop et al. (|2010 | ) see ing a more pronounced bi- modality 
and Chap in et al.l (|201lh a greater tail beyond z = 2. This 
comprises a significantly higher fraction o f high redshift 
sources than in this work. IPve et ahl (|2009l ) find ^ 30% of 
their BLAST sources within a deep field to be at 2; > 1, fully 
consistent with our fraction. Having used different selection 
criteria for their BLAST sources, either signal-to-noise cut- 
offs or flux hmits in any of the 250, 350 or 500 fim BLAST 
bands, leads to different sub-samples that are difficult to 
compare. In addition, the methods to obtain photometric 
redshifts and to identify optical/mid-infrared counterparts 
vary. It is difficult to disentangle the different approaches, 
but there is still very broad agreement in the conclusion 
that we see two different populations, one which, at lower 
redshifts, we find in VIKING counterparts, and the other, 
at higher redshifts, the fraction of which we can imply and 
which is consistent with at least some of the BLAST find- 
ings. 

A similar picture emerges f rom the He r MES p roject 
(lOliver et al.l I2OIOI) so far. "Eal es et al.1 (|2010l i and 
lElbaz et al.l (|2010h use deep imaging in small areas (< 
0.1 deg^) observed by SPIRE and with excellent multi- 
wavelength data available, as well as spectroscopic and pho- 
tometric redshifts. Both groups use a 250 /am selected sam- 
ple and assume a 24 /xm detectio n. With a h i gh fra ction 
of spectroscopic redshifts (> 65%V [Elbaz et all (|2010h find 
35 — 40% of their sources lie at z > 1 (de duced from 
their fig. 2), consistent with our findings, whereas bales et al] 
(|201Q| 1 discover close to 50% in this redshift range. 



5.6 Towards more complete identifications 

So far we have estimated that 73% of SPIRE sources have 
counterparts in VIKING, while 51% have a reliable match; 
thus, the reliable sample comprises approx 51/73 = 70 per- 
cent of all SPIRE galaxies with both /250 > 32 mJy and 
Ks < 19.2. We also know that our blank 14% of SPIRE 
sources are definitely undetected in VIKING, while the 
remaining 35% of SPIRE sources have one or more low- 
reliability match(es); overall we expect around half of these 
to be genuine matches. 

To make this decisive, we would need follow-up observa- 
tions such as radicF^ or ALMA (Atacama Large M illime- 
ter/submillimetre Arrav, IWotten Thompsonll2009l ) imag- 
ing giving a sub-arcsec position, or possibly optical/NIR 
spectroscopy of VIKING candidates (if we can identify 
SPIRE sources via unusually strong emission lines). This 
would normally be "decisive" in that a sub-arcsec ra- 
dio/submm position would either match a VIKING galaxy 



Better positions, and greater efficie ncy of IDs will be p ossible 
with the ASKAP radio survey EMU (|Norris et al.ll201lh which 
will have 10" angular resolution and cover a redshift range quite 
similar to H-ATLAS. 



14 S. Fleuren et al. 



to very high rehabihty, or if not it would prove the source is 
fainter than the VIKING hmit. 

For non-rehab le matches, the total reliability ^ Ri for a 
given source is a good estimate of the probability that the 
SPIRE source has a real counterpart in VIKING, i.e. the 
probability that a follow-up will actually find a good match; 
therefore in a follow-up search we should target the non- 
reliable matches in descending order of total probability. 

Assuming a lower limit of the total probability of 50% 
(70%) would result in an additional 2380 (1856) identifica- 
tions for 2967 (2101) observing targets. We would then ob- 
tain a sample which is 85% (82%) complete to /250 > 32 mJy 
and Ks < 19.2. If we would use all non-reliable positions as 
targets, regardless of the total reliability, we would only be 
able to reach a conipletene ss of 89%. This effect can also be 
seen m ISmith et al.l(|2011al ) where the sum of the reliabilities 
to all possible counterparts would result in a 44% identifica- 
tion rate, lower than expected from their value of Qo = 0.59. 
This shows that our (and their) reliabilities might be under- 
estimated; this is more evidence for a likely non- Gaussian 
positional error distribution which will be addressed in fu- 
ture publications of the catalogue. 

Of additional consideration is our candidate list for multi- 
ple identifications, see section 5.3. They display a total relia- 
bility of 80% by definition and would hence be included in a 
possible target list. It would be interesting to see how many 
of those could be confirmed as true multiple identifications. 

Fig. 1161 shows the distribution of the total reliabilities for 
the non-reliable SPIRE positions together with the number 
of additional identifications we would expect if we followed 
up all SPIRE positions down to 50 or 70 perc ent. 

It is also useful to compare to the results of IPunlop et al.l 
(|2010h : as described above, they identified a much smaller 
sample of 20 BLAST 250 /xm sources, but benefiting from the 
very deep mult i wavelength data in GOODS-S. In their sam- 
ple, all candidate identifications with z < 1.2 are brighter 
than Ks < 21, while all at 2; > 1.2 are fainter; this suggests 
that the 3011 sources without a VIKING counterpart have 
a high probability of being at z > 1.2 , and the same applies 
for the 3194 sources with low-reliability matches ^ R < 0.2. 

Since the flux ratio /sro / /250 strongly increases with red- 
shift for typical SMGs, the non-identified sources are there- 
fore good targets for ALMA 870 /xm follow-up snapshots; 
this could give a relatively efficient method for selecting lu- 
minous high-z SMGs. 



6 COMPARISON TO OPTICAL 
IDENTIFICATIONS 

6.1 Reliable counterparts 

In this section we compare our results with the optical iden- 
tification supplied with the H- ATLAS G09 Phase 1 source 
catalogue. This used a similar likelihood ratio method with 
an r-band select ed sample down to r = 22.4, as explained for 
the SDP data in lSmith et all (l2011al ). The VIKING i^^^-band 
should be better placed than the optical r-band in identify- 
ing counterparts to the SPIRE sources. As discussed in the 
introduction, at higher redshifts (z > 0.5), the Kg band is 
detecting fiux from the near-infrared restframe, while the r 
band is restframe blue/UV; thus Ks is much better able to 
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Figure 16. — Total reliability distribution of the non-reliable 
SPIRE sources. For the 7391 sources without a reliable coun- 
terpart, we sum the reliabilities of all possible VIKING counter- 
parts. This total reliability is an estimate of the probability that 
the SPIRE source has a real counterpart in VIKING. If we include 
all SPIRE positions with total reliabilities greater than 50% for 
a follow-up, we will have to observe 2967 positions and expect to 
gain 2380 additional reliable VIKING counterparts (red line and 
text). The blue line and text show the equivalent for a 70% total 
reliability. 



detect dusty galaxies. We therefore expect a higher num- 
ber of reliable identifications from matching with VIKING 
than with SDSS. To be able to compare our results, we 
cross-match our VIKING candidate matches with the SDSS 
database (DR7) within 2" and choose the nearest (primary) 
object. 

We concentrate on the reliable counterparts of both sur- 
veys. The Phase 1 catalogue lists 9623 reliable optical coun- 
terparts (36.5%) to 26,369 5cr SPIRE sources in the G09 
field; of the reliable counterparts, there are 8587 (39.0%) in 
the VIKING observed area. In comparison, we are able to 
match 11,294 SPIRE sources (51.3%) reliably to VIKING 
Ks objects. 

We find a reliable Ks counterpart for 3732 SPIRE sources 
without a reliable optical counterpart. Of the 3732 positions, 
1717 are blank in SDSS (- 21% of ah SDSS blank fields), 
i.e. they are too faint to be detected in SDSS. Fig. [T71 shows 
the Ks magnitude distribution of the counterparts to the 
1717 SPIRE positions that are optical dropouts (red solid 
line). Unsurprisingly, the magnitudes are rather faint, with 
a median of Ks = 18.26, compared to the magnitudes of all 
reliable Ks counterparts with a median of Ks = 17.07. 

The remaining 2015 SPIRE positions have optical coun- 
terparts, but their reliabilities lie below the threshold of 
R > 0.8. Fig. [17] (black solid line) shows the r modelmag 
distribution of the 3085 candidate matches to those 2015 
SPIRE sources from which we can see that they belong 
mostly to the faint end of the overall magnitude distribu- 
tion. 

Conversely, there are 1,024 sources with optical reliable 
counterparts for which we did not find a reliable Ks counter- 
part. Only 121 of those positions are blank in Ks, mainly due 
to quality issues like saturation or bad pixels; the remain- 
ing 903 sources share 2261 VIKING candidates, of which 
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Table 3. — Comparison of the reliable counterparts to the SPIRE sources using optical r-band and near-infrared Ks-band matching. In 
the area corresponding to the VIKING preliminary source catalogue, we match reliably 51.3% of all SPIRE sources. This is a significant 
increase in the identification rate compared to the 39.0% of sources that are matched reliably to a SDSS object. 
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Figure 17. — The magnitude distributions of reliable Kg coun- 
terparts to SPIRE sources without reliable optical counterparts. 
The black histograms show r (modelmag, AB) distributions, the 
red histograms show Kg (aperture, Vega) magnitude distribu- 
tions. The dashed histograms represent all reliable counterparts. 
The black solid line shows the r magnitude distribution of the 
non-reliable optical matches to SPIRE sources with reliable Kg 
counterpart. The red solid line shows the distribution of Kg mag- 
nitudes of optical drop-outs with reliable Kg counterpart. 
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Figure 18. — VIKING Kg image of HATLAS 
J090550. 5+002216, 30" on the side. The black cross marks 
the SPIRE position. Red circle: reliable SDSS counter- 
part (R = 0.85), black circle: reliable VIKING counterpart 
(R = 0.96). In the optical, the VIKING counterpart is detected 
with r = 23 .16, fainter than th e limit of r = 22.4 used in the 
matching bv lSmith et"all tollaT) . 



706 have reliabilities with 0.5 < R < 0.8. Comparing to our 
multiple candidates, see section [531 we find that 590 of our 
903 sources are indeed included in our candidate list of 1444 
sources. We also find 14 sources that have confirmed multi- 
ple counterparts (by spectroscopic redshift). Over half of the 
reliable matches we miss when compared to the optical iden- 
tifications, could hence be genuine multiple counterparts. 

It is interesting to consider for how many SPIRE sources 
the VIKING and SDSS matching disagree on reliable coun- 
terparts. 7563 SPIRE sources (~ 88% of the reliable optical 
matches in the VIKING area, - 67% of the reliable VIKING 
matches) are matched reliably in both surveys. Here, 7404 
are matched to the same object. This leaves only 159 SPIRE 
sources (2.1% of matches) where the identification disagrees. 
Some of those are deblending issues; often though we find 
that the reliable optical counterpart is too faint in the Kg 
band and/or the VIKING counterpart is too faint or not 
detected in the r-band, resulting in different identifications. 
Fig. [m shows an example, HATLAS J090550. 5+002216. 

6.2 Stellar counterparts 

Of our 12 reliable stellar matches, 7 have reliable SDSS 
counterparts. Here, 2 have been classified as galaxies in the 
Phase 1 catalogue, HATLAS J091233.9-004549 and HAT- 



LAS J085353. 2+001648 (10930 and 21662 in Fig. Both 
are clearly stellar, showing diffraction spikes in both SDSS 
and VIKING images. The reliabilities of the stellar matches 
to the remaining 5 sources are high with 0.43 < R < 0.78. 

Conversely, in the optical catalogue there are 21 reliable 
stellar matches: of these, 19 are in the VIKING area, of 
which 5 are matched reliably to a VIKING star, 7 stellar 
objects are not included in our VIKING sample due to sat- 
uration in Kg, and 7 more are matched, but do not reach 
R > 0.8 (but have reliabilities > 0.4). 

6.3 Photometric redshift comparison 

Due to the facts that we match with SDSS to obtain our 
photometry and that ANNz is used in both cases, it is no 
surprise that the photometric redshift distributions of the 
reliable counterparts are broadly similar, see Fig. 1191 We 
obtain a slightly higher median redshift of 5 = 0.396 com- 
pared to the median of z = 0.326 of the photometric red- 
shifts supplied with the Phase 1 catalogue, partly due to 
a higher number of redshifts Zphot > 1. Of the 309 SPIRE 
sources with VIKING reliable counterparts and Zphot > 1, 
76 are reliable, have photometric redshifts and are matched 
to the same object in the optical catalogue. The photomet- 
ric redshifts of those 76 objects differ by an average of 0.46. 



16 S. Fleuren et al. 



2500 



2000 



1500 



1000 



500 



^Dhot' Phase 1 catalogue 
Zphot,all VIKING 
Zphot, VIKING with SDSS 



0.0 



0.4 



0.8 



1.2 



1.6 



Figure 19. — Comparison of the photometric redshift distri- 
butions obtained in this work (VIKING) and supphed with the 
Phase 1 catalogue. The red hue shows the latter, 9583 photomet- 
ric redshifts with a median oi z = 0.326. The black dashed line 
shows the former, 11,294 photometric redshifts with a median of 
z = 0.396. The black line represents photometric redshifts of the 
8748 VIKING objects with SDSS counterparts. 



ric redshifts from SDSS/UKIDSS LAS photometry for the 
HATLAS SDP field. For z > 0.3 we find that photometric 
redshifts are systematically underestimated with a bias of 
- 0.09. 

Comparing our results with that from the r-band match- 
ing to SDSS objects supplied with the SPIRE catalogue, we 
report a > 12% increase in the reliable identification rate. 
We find that we agree on reliable counterparts for ^ 88% of 
the reliable optical matches to sources within the VIKING 
area of the G09 field. 

The identifications here provide a useful potential pre- 
selection for follow-up studies: the moderate-reliability 
matches could mostly be confirmed or rejected using optical 
multi-object spectroscopy, giving a mostly complete subsam- 
ple for sources at z < 1; while the ^ 28% of sources with no 
match or low-reliability match (es) have a high probability 
of being at 2; > 1, and form a large sample of interesting 
targets for ALMA 870 /xm snapshots. 

A future SPIRE source catalogue will include the VIKING 
ZYJ}iKs photometry (aperture magnitudes) for all candi- 
date matches, but we stress that only matches with R > 0.8 
should be regarded as reliable counterparts and used for sci- 
ence application. 



This large difference could be explained by incomplete pho- 
tometric information from UKIDSS LAS used to compile 
the photometric redshifts in the optical catalogue. Indeed, 
nearly half (35) of the 76 objects have only 1 or 2 bands 
in the near-infrared available from LAS. This shows the ad- 
vantage of the deeper VIKING data. Much better results 
should be achieved once we have optical photometry from 
the VST KIDS survey to combine with VIKING. 



7 CONCLUSIONS 

We have matched the 22,000 SPIRE G09 sources that fall 
within the VIKING observed area to a catalogue of near- 
infrared objects from VIKING: the VIKING sample con- 
tains 1,376,606 objects, classified into 847,530 galaxies and 
529,076 stars according to shape and c olour parame t ers us - 
ins; a modified version of the method of iBaldrv et"all (l2010h . 
We found statistically, using blank- field comparisons, that 
73 =b 2% of SPIRE sources should have one or more counter- 
part detections in VIKING. With a search radius of 10^^ we 
found 35,800 candidate matches: applying a likelihood ratio 
method to calculate the probability of each candidate to be 
the true counterpart of the SPIRE source, we find matches 
to 11,294 sources (51.3%) with a probability > 80%, or reli- 
ability R > 0.8. The false identification rate is estimated to 
be 4.2% and the probability of mis-identifying a true coun- 
terpart is ~ 5%. 

Of the reliable counterparts, 3147 (27.9%) have spectro- 
scopic redshifts from either the GAMA or the SDSS redshift 
surveys. We calculate photometric redshifts for the remain- 
ing possible matches, using a sample of 32,465 spectroscopic 
redshifts as a training set. The errors in the redshift es- 
timation are investigated using the existing 3147 spectro- 
scopic redshifts. We find a scatter of cr = 0.0353 in the dif- 
ference \z svec — ^p/^otl which I S Comparable with a = 0.037 
found by ISmith et al.l (l2011al ) when calculating photomet- 



8 ACKNOWLEDGEMENTS 

Herschel is an ESA space observatory with science instru- 
ments provided by European-led Principal Investigator con- 
sortia with significant participation from NASA. U.S. par- 
ticipants in Herschel ATLAS acknowledge support provided 
by NASA through a contract issued from JPL. 

GAMA is a joint European- Australasian project based 
around a spectroscopic campaign using the Anglo- 
Australian Telescope. The GAMA input catalogue is based 
on data taken from the Sloan Digital Sky Survey and the 
UKIRT Infrared Deep Sky Survey. Complementary imag- 
ing of the GAMA regions is being obtained by a number of 
independent survey programs including GALEXMIS, VST 
KIDS, VISTA VIKING,WISE, Herschel- ATLAS, GMRT 
and ASKAP providing UV to radio coverage. The GAMA 
website is: http://www.gama-survey.org/ This work used 
data from the UKIDSS DR5 and the SDSS DR7. The 
UKIDSS project is defined in Lawrence et al. (2007) and 
uses the UKIRTWide Field Camera (WFCAM; Casali et al. 
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APPENDIX A: VIKING CUT-OUTS 

We show a selection of cut-outs around VIKING reliable counterparts to SPIRE sources, drawn at random from 16 magnitude 
bins ranging from Ks magnitude (Vega) 13.4 — 19.8. The cut-outs are 1 arcmin on the side. 



